Hydrodynamically-driven colloidal assembly in the thin-film entrainment regime 
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We study numerically the hydrodynamics of dip coating from a suspension and report a mechanism 
for colloidal assembly and pattern formation on smooth and uniform substrates. Below a critical 
withdrawal speed of the substrate, capillary forces required to deform the meniscus prevent colloidal 
particles from entering the coating film. Capillary forces are overcome by hydrodynamic drag only 
after a minimum number of particles organize in a close-packed formation within the meniscus. Once 
within the film, the formed assembly moves at nearly the withdrawal speed and rapidly separates 
from the next assembly. The interplay between hydrodynamic and capillary forces can thus produce 
periodic and regular structures within the curved meniscus that extends below the withdrawn film. 
The hydrodynamically-driven assembly documented here is consistent with stripe pattern formations 
observed experimentally in the so-called thin- film entrainment regime. 



PACS numbers: 47.85.-g; 82.70.Dd; 47.61.Jd 

The ability to build crystalline microstructures on 
smooth substrates enables significant advances in mate- 
rials science and microfabrication, e.g. optoelectronics 
and metamaterials [TJ |2]. Film deposition from a col- 
loidal suspension, via substrate withdrawal and/or sol- 
vent evaporation, can produce patterned microstructures 
that extend over millimeter lengths [3H6]. For example, 
experimental studies have demonstrated a critical with- 
drawal speed below which the periodicity and regularity 
of the structures increase dramatically [5j [6] . The criti- 
cal speed corresponded closely to the conditions at which 
the particle diameter is equal to the thickness of a thin 
coating film predicted by the Landau-Levich theory [7]. 

These conditions set the upper limit of the so-called 
thin- film entrainment regime. In this regime, the for- 
mation of periodic stripe patterns has been observed on 
partially wetting substrates [5] and further studied on hy- 
drophilic substrates [6]. These experimental studies re- 
port that (i) very low particle concentrations are required 
for the formation of regular assemblies, and (ii) varying 
the withdrawal speed has a much greater effect on the 
stripe morphology than varying the solvent evaporation 
rate. The employed technique, referred to as convec- 
tive self-assembly (CSA), involves convective transport 
of colloidal particles from the bulk of the solution to the 
meniscus and coating film. In regions where the assem- 
bly process takes place, the energies of convective mo- 
tions and capillary interactions are much larger than the 
thermal energy /c#T (here ks is Boltzmann's constant, 
and T is temperature) and Brownian effects may be ne- 
glected. Self-organization processes in CSA can thus in- 
volve nontrivial hydrodynamic interactions between col- 
loidal particles. Previous studies, however, attributed 
the stripe pattern formation to different mechanisms that 
only consider capillary interactions and quasi-static pro- 



cesses [5j[6]. In this Letter, we study and document the 
critical role of hydrodynamic effects in the particle as- 
sembly and subsequent pattern formation. 

To study the hydrodynamics of the fluid-particle sys- 
tem, we employ a class of lattice Boltzmann (LB) meth- 
ods for isothermal multiphase flow [8] [9]. The solid 
phase, which includes a moving wall (i.e. the with- 
drawn substrate), three stationary walls, and N p col- 
loidal particles, is modeled with an immersed boundary 
(IB) approach [10]. The employed IB approach repre- 
sents the solid phase with a smooth distribution function, 
$s(x,£) G [0,1] [see Fig. 1(a)] that determines local hy- 
drodynamic forces due to (short-range) fluid-solid inter- 
actions (e.g. molecular collisions). Long-range molecular 
interactions that give rise to capillary forces are mod- 
eled by two mean- field potentials, ?/; FF (x, t) for fluid- fluid 
and iJj FS (x,t) for fluid-solid interactions [TQ [12]. The 
interaction potential ip FF models a volatile fluid that 
separates into two phases, which are modeled with an 
ideal equation of state [12]. The coexistence densities 
are pL = m/(Ax) 3 and py = 0.1 m/(Ax) 3 for the liquid 
and vapor phase respectively; hereafter, m is the molec- 
ular mass and Ax is the numerical grid spacing. The 
resulting surface tension is 7 ^ 0.13 ksT / (Ax) 2 . We ad- 
just fluid- solid interactions so that the solid phase is al- 
ways wetted by the liquid phase, which prevents colloidal 
particles from breaching the liquid- vapor interface. For 
this work, we neglect long-range molecular interactions 
between solid bodies. We also neglect thermal fluctua- 
tions that produce Brownian effects. The modeled fluid- 
solid forces, F FS = T [$5, ip FF ip FS ] , govern the transla- 
tional and rotational dynamics of the colloidal particles. 
The equations of motion for the center-of-mass position, 
x p (t), and angular velocity, u* p (t), of each particle p are 
integrated using a conventional leapfrog scheme. The nu- 
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merical method is described in detail in the supplemental 
material. 

The numerical model is applied to study the dynamics 
of colloidal assembly in the thin- film entrainment regime. 
Dip coating of plates produces a nearly two-dimensional 
flow that drives the convective transport of particles from 
the suspension bulk to the meniscus region and subse- 
quently to the substrate. Furthermore, the formation 
of regular stripes [2 EJ [6] requires a statistically homo- 
geneous distribution of particles along the direction of 
the stripes. We proceed by studying a two-dimensional 
flow configuration that resembles the classical Landau- 
Levich problem [7] of plate coating via withdrawal from 
a liquid bath, where we later introduce circular parti- 
cles. Notwithstanding these geometrical simplifications, 
we capture the basic mechanisms of particle assembly 
and stripe formation. 

We perform numerical simulations in a rectangular do- 
main that is partially bounded by walls. The domain size 
is L x x L y ~ 7 £ c x 3.2 £ c where the capillary length 
is £ c — \Zj/(pl ~ Pv)9 — 133 Ax with a gravitational 
acceleration magnitude g = 5 xlO~ 6 Ax / (At) 2 . At ini- 
tialization, the liquid bath lies below yn — 2.25 £ c (see 
Fig. 1). The vertical speed on the left wall, i.e. the with- 
drawal speed, is prescribed to small values, U ~ 0.25- 
2.0 xlO _2 Ax/At, which results in moderately low capil- 
lary numbers Ca = Up/j ~ 0.01-0.08. Neutrally buoy- 
ant particles of circular shape and radius R = 10 Ax are 
used for the simulations in this Letter (the supplemen- 
tal material includes results demonstrating grid conver- 
gence). The particle dynamics are thus characterized by 
small Reynolds numbers Re = piJJRj p ~ 0.075-0.3 and 
a small Bond number Bo = (R/£ c ) 2 = 6 xlO -3 . Parti- 
cles in varying numbers, from zero to twenty, are released 
at the same time to with zero initial velocity from differ- 
ent locations within the left half of the liquid bath; i.e. 
x p (t ) e (x w +R,L x /2] and y p (t ) G [vh -20R,y H -3R] 
(xw indicates the surface of the coated plate). Details of 
the simulation procedure are provided in the supplemen- 
tal material. 

First, we study the background flow topology. Simu- 
lations without particles reach steady state after a time 
ts = (L y — yH)/U. The steady-state flow topology, re- 
ported in Fig. 1, is typically observed in dip coating [15] . 
A fundamental feature is the presence of a stagnation 
point S* on the vapor-liquid interface at vertical posi- 
tion y = y* [cf. Fig. 1(a)]. The flow streamline ending 
at S* defines the boundary between a "shear flow" re- 
gion, Tl, which continues into the coating film, and a 
"recirculation flow" region, T2, which extends into the 
liquid bath. Under the studied conditions, the entrained 
film is always thinner than the particle diameter; one 
simulated particle is shown in Fig. 1(a) for reference. 
The thickness profile h(y) of the dynamic meniscus be- 
low the film largely determines whether a particle can 
enter the film. Above the meniscus (for y > yn + 2£ c ), 
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FIG. 1. Flow topology. (a) Streamlines in the menis- 
cus (Ca=0.04). A stagnation point S* (circle) lies on the 
interface where h(y*)=h* . The streamline (dashed line) 
ending at S* separates a "shear flow" region, IT, from a 
"recirculation flow" region, T2. The surface of the with- 
drawn substrate is at x—xw and the horizontal level of 
the bath (at initialization) is at y—yH- A particle of size 
i?=10Ax is illustrated for reference; the diffuse fluid-solid in- 
terface is centered at $s(x, t)=0.5. (b) Normalized thickness 
h* /2R=(h* li c ) I (2VBo) vs. Ca=Ufi/j. Simulation results 
(symbols) compare well against lubrication theory (solid line) 
[7]. The critical entrainment condition h*—2R corresponds to 
Ca*~0.014. 

the film thickness approximates the analytical prediction, 
hf/l c = 0.95 Ca 2 / 3 - 0.1 Ca, valid for Ca < 0.1 [7 . More 
importantly, the meniscus thickness h(y*) = h* observed 
at the stagnation point agrees well with predictions based 
on lubrication theory [7]. In Fig. 1(b) we report this char- 
acteristic thickness normalized by the particle diameter 
h*/2R = [3(h f /£ c ) - (h f /£ c ) 3 /Ca]/2VBo~. There is a 
withdrawal speed corresponding to Ca = Ca* = 0.014 
for which h* = 2R [see Fig. 1(b)]; this condition appears 
to be critical for particle entrainment as we discuss below. 

We now proceed to analyze the trajectories 
{x p (t), y p (t)} after the background flow is established. 
For this purpose, we release particles simultaneously 
at to = 1.2 ts- We begin with the simplest case of 
a single particle, which yields a variety of open and 
closed trajectories. Single particles released within 
the "recirculation flow" region T2 attain stable orbits 
within the bath that prevent them from entering the 
film (see supplemental material). The case of a single 
particle (p=l) released within the "shear flow" region 
Tl is of particular interest. When approaching the 
film, the trajectories of single particles within Tl are 
compressed into a unique "entrainment" trajectory 
by volume exclusion effects. Consequently, when the 
withdrawal speed is sufficiently large, i.e. for Ca > Ca* , 
a single particle successfully enters the film. The 
vertical displacements Ay = yi(t) — yu > for different 
withdrawal speeds corresponding to Ca = 0.02-0.08 
are reported in Fig. 2(a). A significant variation in the 
particle speed y = dy/dt can be observed in the inset of 
Fig. 2(a): the speed is nearly constant before entering 
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FIG. 2. Single particle dynamics for Ca > Ca*. (a) Verti- 
cal displacement Ay — yi(t) — yH vs. tU/R for Ca > Ca* . 
The time origin (t = 0) is chosen so that Ay(0) = 0. For 
Ay > individual particles released within Tl follow a com- 
mon xy-trajectory at a given Ca. The inset in (a) shows the 
vertical speed yi — dyi/dt. (b-e) Streamlines for Ca = 0.04 
at different time instances indicated in panel (a). 



the meniscus (yi/U c± 0.65 for Ay/R < 10); decreases 
inside the meniscus (10 < Ay/R < 20); and then 
increases after clearing the meniscus and entering the 
flat film (Ay/R > 20). The minimum vertical speed is 
reached near the stagnation point that is initially located 
at y* . While the minimum speed is directly correlated 
with the confinement ratio h* /2R, the maximum vertical 
speed reached within the film is inversely correlated 
with h*/2R. The vertical speed is dominant along the 
entrainment trajectories in Fig. 2(a); i.e. \xi\/\yi\ < 0.1 
and \wiR\/\yi\ < 0.75 (wj is the angular velocity). 

We next focus on the force balance along the nearly 
one-dimensional trajectory for entrainment. As seen in 
Figs. 2 (b-e), the entrainment of a single particle is ac- 
companied by significant deformation of the meniscus 
and the flow streamlines. In the horizontal direction, 
capillary forces are mainly balanced by a high pressure 
within the particle- wall gap and xj <^U. In the vertical 
direction, hydrodynamic drag driving the particle motion 
is counteracted by capillary forces. For a confined par- 
ticle moving parallel to the wall, the Stokes drag may 
be expressed as Fp = —fi(Kiyi — KuU), where Kj > 
and Kjj > are resistance coefficients. The ratio of re- 
sistances fr> = Kjj/Kj — >• for low confinement and 
far from the wall (h/2R ^> 1), while fp ~ 1 in high 
confinement conditions and close to the wall (h ~ 2R). 
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FIG. 3. Single particle dynamics for Ca < Ca* . (a) Vertical 
displacement Ay = yi(t) — yH vs. tU/R. The time origin 
(t = 0) is chosen so that Ay(0) = 0. Single particles re- 
leased within n reach a steady position inside the meniscus 
for tU/R>25. (b-c) Streamlines for Ca = 0.01 < Ca* at two 
time instances indicated in panel (a). 



The capillary force (per unit length) is Fq = —jdA/dyi 
where A is the total interfacial area (per unit length). 
For Bo = (R/£c) 2 <C 1, we consider that capillary forces 
act normal to the interface with their magnitude deter- 
mined by a shape function fc(h/2R) = \dA/dyi\. We 
then neglect inert ial effects of order 0(Re) to obtain 
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U 



f D (h/2R) 



fc(h/2R) dh/dy 



KjCa + \dh/dy\* 



(1) 



where fjj and fc are positive functions determined by 
the thickness profile h(y). Hence, Eq. ([!]) predicts a decay 
in vertical speed within the meniscus where dh/dy < 0; 
this capillary effect increases for small Ca. Once inside 
the film, we have dh/dy ~ while the high confinement 
increases the drag, fo(h/2R) — ^ 1, and thus the particle 
speed yi — » U. 

Moreover, the force balance in Eq. predicts that 
colloidal particles can be trapped (yj = 0) within the dy- 
namic meniscus below a critical capillary number Ca*. 
Indeed, simulation results in Fig. 3 for a low withdrawal 
speed, where Ca < Ca* = 0.014 and h* < 2R, show that 
a single particle attains a stationary position within the 
meniscus. Thus, Eq. indicates, through the functions 
fn and fc, that widening and/or flattening the menis- 
cus has critical effects on the particle speed. For single 
particles we observe that yj — >• for h* /2R — > 0, while 
yi ~ const, for h* /2R — > oo (cf. Figs. 2-3). This observa- 
tion suggests that a leading particle widens the meniscus, 
which would allow particles behind to climb faster. This 
scenario is graphically illustrated by the slope and in- 
tercept of the straight lines in Fig. 2(a) for h* /2R = 1.4; 
particles within 9R below a leading particle would assem- 
ble in the meniscus if moving with the speed (y 0.65 U) 
observed for h* > 2R. Moreover, particles already in the 
film would climb away about 26.5 R during the assembly. 

Simulations with multiple particles confirm a 
hydrodynamically-driven assembly within the meniscus 
for Ca = 0.01 < Ca*. Three particles (p=I-II-III) in one 
case [cf. Fig. 4(a)] and two particles (p=I-III) in a second 



4 




FIG. 4. Hydrodynamically-driven assembly for Ca = 0.01 < 
Ca*. (a-b) Vertical displacement Ay p = y P (t) — y P (0) vs. 
tU/R. Three (p=I-II-II) and two (p=I-II) particles are re- 
leased at t = 0; yjj(0) = yi(0) - 3R, ym(0) = yi(0) - 6R, 
and yi(0) ~ (c-e) and (f-g) Streamlines at time instances 
indicated in panels (a) and (b). 



FIG. 5. Multiple particle assembly in thin films. A sequence 
shows six time instances for (a) Ca=0.01 and 0=0. 13R, (b) 
Ca=0.01 and 0=0. 25R, and (c) Ca=0.04 and 0=0. 13R. For 
Ca < Ca* and low area factions < 0.13, the assembly occurs 
in the dynamic meniscus, which increases the periodicity and 
regularity of the arrays (see movies in supplemental material) . 



case [cf. Fig. 4(b)] are released simultaneously with 
initial separations yi — yn = 3R and yi — ym = 9R 
where yi(0) c± y^. The vertical displacements in 
Fig. 4(a) show that the trailing particle (III) climbs 
faster than the leading particles (I-II) , as hypothesized 
above. The intercept of the straight lines in Fig. 4(a) 
estimates maximum initial separations of 6-8R for which 
the assembly is possible. In both cases the particles 
reduce their initial separations and assemble inside the 
meniscus [cf. Figs. 4(c-e) and Figs. 4(f-g)]. Below 
the critical condition Ca < Ca*, a minimum number 
of colloidal particles (three in the studied conditions) 
must assemble in a close-packed formation before film 
entrainment is possible. A close packing of particles 
increases the hydrodynamic drag and flattens the local 
interface curvature so that dh/dy —> for particles in- 
side the assembly. When hydrodynamic drag overcomes 
interfacial forces, the particle assembly enters the film 
and increases its speed as seen in Figs. 4(c-e). 

Simulations with larger numbers of particles (N p =10- 
20) further corroborate the periodic formation of reg- 
ular assemblies when Ca < Ca* [cf. Fig. 5(a)]. To 
draw this conclusion, we performed five realizations of 
each numerical experiment by releasing particles from 
random locations with root-mean-square particle sepa- 
rations |Ax rms | = 3.5R and |Ax rms | = 2.5i?, which cor- 
responds to area fractions = 0.13 and = 0.25 re- 



spectively. The observed assembly within the dynamic 
meniscus was periodic and nearly identical (see movies 
included with the supplemental material). For Ca < Ca* 
and the lower area fraction = 0.13 [cf. Fig. 5(a)], par- 
ticle assemblies fully form before entering the film, which 
increases the periodicity and regularity of the arrays. As 
seen in Figs. 5(b-c), increasing the area fraction or the 
withdrawal speed produces irregular assemblies that form 
once within the coating film. In the latter case, a differ- 
ent assembly mechanism, due to long-ranged capillary 
interactions and film instabilities PQ, deteriorates the pe- 
riodicity and regularity of the observed stripe patterns. 

In conclusion, our study indicates that below a 
critical withdrawal speed, for which h* < 2R, a 
hydrodynamically-driven assembly within the meniscus 
can produce highly organized structures. The identi- 
fied mechanism considers the coupling between hydro- 
dynamic and capillary interactions. Our simulations and 
analyses qualitatively agree with recent experimental ob- 
servations of the formation of regular stripe patterns 
[5j [6] . The results in this Letter lead to a better funda- 
mental understanding that could improve current tech- 
nologies used for the flow-driven assembly of colloidal 
crystals. 
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